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Abstract. Convective instability in a ferrofiuid layer heated from below or from above 
in the presence of a uniform vertical magnetic field is investigated theoretically. Con- 
vection is caused by a magnetic mechanism based on the temperature and concen- 
tration dependence of magnetization. An imposed temperature gradient establishes 
(by the Soret effect) a concentration gradient of magnetic particles of which the fer- 
rofiuid is composed. Both these gradients cause a spatial variation in magnetization, 
which induces a gradient of magnetic field intensity within the fluid layer. The field 
gradient induces in its turn an additional redistribution of magnetic grains due to mag- 
netophoresis. Resulting self-consistent magnetic force tries to mix the fluid. A linear 
stability analysis predicts oscillatory instability in a certain region of the magnetic field 
strength and the fluid parameters. The instability owes magneto- and thermophoresis 
its origin: were the particle diffusion not operative, then only stationary instability 
would occur. A discovery of predicted convective oscillations is expected in ferrofiuid 
layers about 1 mm thick, where the buoyancy mechanism is negligible and the charac- 
teristic diffusion time is not too long. 

1 Introduction 

Magnetic colloids (so-called ferrofluids [1,2]) represent very interesting object for 
studies of convection. Even in the absence of the gravity, mechanical equilibrium 
of a nonisothermal ferrofiuid in a magnetic field H is in general impossible. At the 
basis of the mechanism of thermomagnetic convection [1-5] lies the temperature 
dependence of the magnetization M: under otherwise equal conditions, the colder 
elements of the fluid are more strongly magnetized, and therefore they are subject 
to a larger magnetic force in the direction of VH. The gradients of the magnetic 
field intensity here play the same role as does the gravitational field g in the 
mechanism of ordinary thermogravitational convection, whilst the temperature 
dependence of the magnetization, M(T), plays the role of T-dependence of the 
fluid density p. 

Interestingly, thermomagnetic convection can arise even in a uniform applied 
magnetic field [4] . Let a ferrofiuid confined between two horizontal planes be in a 
uniform vertical external field H e = (0, 0, H e ) at the temperature T(z). Then the 
dependence M(T) leads to the result that the field inside the fluid, H = (0, 0, H), 
appears to be nonuniform. Indeed, from the Maxwell equation div B = and 
the definition of magnetic induction B = H + 47rM it follows 
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that yields 



dH 

d7 




(1.2) 



where fi = 1 + 4n(dM/dH) is the differential magnetic permeability. Equation 

(1.1) shows that the induced gradient of magnetic field is directed always opposite 
to the gradient of magnetization. It results in unstable magnetization stratifi- 
cation: the magnetic force F = MVif tries to mix ferrofluid. According to eq. 

(1.2) , this force is proportional to the imposed temperature gradient. Therefore, 
when both the magnetic and buoyancy mechanisms are operative, the first of 
two predominates over the second one at the onset of convection in a thin enough 
fluid layers, where the critical temperature gradient is sufficiently large. One can 
show [1] that the condition 



which permits neglect the buoyancy mechanism of convection, is reduced to 



where d is the layer thickness, Ra c ~ 10 3 the critical Rayleigh number (Ra c = 
1708 in the case of rigid boundaries and Ra c = 277r 4 /4 w 657.5 for free bound- 
aries), n and k the fluid viscosity and thermodiffusivity. For commonly used 
magnetite ferrofluids on the base of water or kerosene the pyromagnetic coeffi- 
cient -(dM/dT)u in the field H ~ 200 - 500 Oe is about 0.50 G/K where 4> 
is the volume fraction of magnetic grains. Then for ~ 0.1, r\ ~ 10 _2 Ps, k <~ 
10- 3 cm 2 /s, and -(dp/dT) p - 5 x lO^gcm^K" 1 , the inequality (1.3) holds 
true down to d = 1 mm. The same conclusion has been done by Finlayson [4] 
who first studied theoretically convective instability of a ferrofluid layer heated 
from below in the presence of a uniform vertical magnetic field. 

It is important to notice that in [4] and posterior theoretical works [5-8] fer- 
rofluids were treated as pure (i.e., single-component) magnetized fluids. Such 
a model assumes a fixed uniform concentration of magnetic grains, <j> — const. 
Ignoring in that way diffusion effects, authors of [4-8] came to conclusion that 
only stationary convective instability is possible in ferrofluids. As it pointed 
out already in the abstract of Ref. [4], "oscillatory instability cannot occur". 
But a real ferrofluid should be treated as a binary mixture with allowance for 
concentration dependence of its magnetization. An externally imposed temper- 
ature gradient induces in a ferrofluid layer a concentration gradient owing to 
the Soret effect. Both these gradients cause a spatial variation in magnetization 
M{<p, T, H) that leads in turn to appearance of a magnetic field gradient within 
the fluid layer. Thus there arise magnetic forces, which try to mix the ferrofluid 
and, in addition, they induce a certain redistribution of magnetic grains due to 
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(dM/dT) H 

9(dp/dT) p 



(1.3) 
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magnetophoresis. As it will be shown below, just the magnetophoresis together 
with the thermophorcsis (the Soret effect) give rise to oscillatory instability in a 
certain region of the magnetic field intensity and the fluid parameters. Appro- 
priate experimental conditions to observe predicted instability take place in a 
thin fcrrofluid layer, when the ordinary (buoyancy) mechanism of convection is 
negligible and the characteristic diffusion time, td = cP/n 2 D, required to settle 
a steady state, is not too long (here D is the diffusion coefficient of magnetic par- 
ticles in the basic liquid). Recent measurements [9] yielded D = 2.7 x 10~ 7 cm 2 /s 
in water and cyclohexane. (This magnitude of D agrees well with a theoretical 
estimate by the Einstein's formula for the diffusion coefficient, D — knT/Qirqa, 
where 2a ~ 10 nm is the diameter of a magnetic grain). According to above for- 
mula for td, the concentration equilibrium in a fcrrofluid layer of the thickness 
d = 1 mm is reached after one hour. 

In the next Section we obtain the description of the magnetophoresis and 
Soret effect in a fcrrofluid assuming that the fluid behaves with respect to an 
external magnetic field like a Langevin's paramagnet [1,10]. In Sees. 3 and 4 we 
derive the complete set of modified linear equations for perturbations and settle 
the boundary conditions for them. The eigenvalue problem is solved in Sec. 5 
for the Benard configuration. On the base of the solution, oscillatory instability 
is predicted to occur for a wide range of magnitudes of magnetic field, Soret 
coefficient, and other fluid parameters. 

Owing to the small magnitude of D (see above), there are possible two sce- 
narios of the beginnings of convection. Depending on the imposed heating rate, 
the fcrrofluid may behave cither as a binary mixture or like a pure fluid. Both 
the possibilities are discussed in Sec. 6. 



2 Magnetophoresis and Soret Effect 

In a nonuniform magnetic field there acts on each particle a force f — (mV)H. 
Orientation of the particle magnetic moment m in the direction of H is impeded 
by thermal motion. As a result, the mean value of m appears to be equal to 

(m) = m£ (S) e ' e = § 5 ^) = cothe-r 1 - (2.i) 

With allowance for the identity H x rotH = ^Vi? 2 — (HV)H and the condition 
of absence of electric currents in nonconducting ferrofluids, rotH = 0, formula 
(2.1) enable us to write the mean magnetophoretic force (f) = m£(£)Vi7. This 
force and the Stokes drag coefficient, 6nr]a, for a sphere of the radius a determine 
the regular component of the Brownian velocity of the particle with respect to 
the liquid, 

U = ^AOVH. (2.2) 

The volume density of the diffusion flux of the matter in a ferrofluid may be 
written 

j = 0u-^(V(/. + 5tVT), (2.3) 
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where St is the Soret coefficient. Substituting (2.2) in (2.3) and using the Ein- 
stein's formula for the diffusion coefficient, we obtain 



j = -D[V(t> + StVT - (<f>/H)t£(OVH}. 



(2.4) 



Thus the matter flux is provided by the " diffusiophoresis" (oc V</>; this term was 
proposed by Derjaguin et al. [11]), thermophoresis (oc VT), and magnetophoresis 
(oc VH). If, however, confined horizontal planes are impervious and an applied 
magnetic field is uniform, then gradients of and H can be caused only by an 
imposed temperature gradient. Therefore in the steady state (j = 0) the values 
of V(j) and VH are expressed via VT. As it is seen from (2.4), the equilibrium 
equation j z = reads 



dT 

'cb 



ding 

dz 



(2.5) 



Substituting in (1.1) M — cf)M s C{mH / k^T) where M s is the saturation magne- 
tization of the particle material, we find 



dln<£ £C'dlnT 



dz 



C dz 



(2.6) 



Eliminating now the concentration gradient from (2.5) and (2.6), we arrive at 

AH 4nM 



dz 



aT 




(2.7) 



where a = fi + UttxqC 2 ^) = 1 + 12n X o(C' + £ 2 ), x = <j)M s m/3k B T the 
initial Langevin magnetic susceptibility, and \P = {T /4>)St plays the role of the 
separation ratio in the mechanism of thermomagnetic convection. Substituting 
dH/dz from (2.7) into (2.5) we obtain, finally, 



cb 



P4_ 
aT 



(A-i) 



dT 

d7 



One can rewrite the latter in the form 

dz T H dz' 
introducing the field-dependent separation ratio 



(A-i) 



(2.8) 



(2.9) 



(2.10) 



Note that <Ph differs from zero even if & = 0, i.e., in the absence of the Soret 
effect. The non-Soret <j) — T coupling originates from the magnetophoresis, which 
is manifested as a negative separation ratio: 



^ H k=o = -(£-i)£C/<7<o. 



(2.11) 
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Actually, the induced gradient of magnetic field (2.7) is directed along the tem- 
perature gradient, thus the hotter fluid (and hence the stronger magnetic field) 
- the higher concentration of magnetic grains. By the same manner it behaves 
a nonmagnetic binary mixture with a negative separation ratio. The function 
defined by Eq. (2.11) is proportional to £ 2 at £ <C 1 and inversely propor- 
tional to £ at £ 3> 1. For a ferrofluid with the initial magnetic permeability 
^0 = 1 + 47rxo = 3 this function has a distinct minimum, <Ph|#=o = —0.303, at 
£ = 2.30. 



3 Equations for Stability Analysis 

The general expression for the volume density of magnetic force acting a non- 
conducting and incompressible ferrofluid may be written [12] 

+ MVH, (3.1) 





where \ = M/H stands for the magnetic susceptibility. The first term in (3.1) is 
known as the magnetostrictive force and the second one called the Kelvin force. 
Under the substitution F in the equation of fluid motion, the magnetostrictive 
force which originates from a potential is included in the pressure gradient, 
Vp, so this force is equilibrated automatically. Thus it remains only the Kelvin 
force, which in equilibrium should be also balanced at each point by the pressure 
gradient: Vp = MVH. On applying the operation rot to this equation, we obtain 
a necessary condition for equilibrium, 

In the case under consideration, all three gradients contained in (3.2) are parallel 
to each other (they have only z-components), so the condition (3.2) is satisfied. 
Thus, the equilibrium is possible, and we set about a study of its stability. 

Small perturbations of a standing mode will be characterized by velocity 
v, pressure p, temperature 9, concentration <p, and magnetic field h. Then the 
perturbation F' of the magnetic force density F = MVH is expressed via 9, 
and h: 

^-^M + (f, + ^ + ^)v ff . (3.3, 

[We have omitted here the term V(Mh z ). Later on it is assumed to be included in 
the gradient of the pressure perturbation, Vp, see Eq. (3.6)]. In the equilibrium 
VM and ViJ are directed along the z axis; so, using (1.1) we obtain from (3.3) 
the only component of the vector F': 
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Substituting in the formula M — 0M S £(£), £ = mHjk^T, and dH/dz from 
(2.7), we get 



/iftz + 47rM 



^ £ T) 



dT 

dz 



Thus, the linear equation of fluid motion takes the form 



dv 
"dt 



p— = -Vp + ??V 2 v+ , | W + 



MAT ( , ££' 



47rM 



(3.5) 



e. (3.6) 



We have substituted in (3.6) the magnitude of the stationary temperature gradi- 
ent, dT/dz = AT Id, where AT is the temperature difference between confined 
surfaces of the layer. In (3.6) e is the unit vector along the z axis, and <P the 
potential of magnetic field perturbations, h = V<£. To derive an equation for 
one should substitute the perturbation of magnetic induction, 



M 

b = h + 4?r— (h - h z e) + 4ir \ 



dM „ dM 



dM . 



in the Maxwell equation divb = 0. It yields 



, ^d 2 <P , / 1 dip ££' ae\ 

^-^-^ = -^ M {^-TXd- Z )- 



(3.7) 



where fi=l + 4ir{M/H) = 1 + 127rxo£(£)/£- 

Equation (2.4) under the replacement <f> — > 4> + <p, T — > T + 9, and H — ► 
i? + d<P/dz determines the matter flux perturbation. So, making allowance for 
Eq. (2.9), we obtain the diffusion equation 



dp 
~di 



Equations (3.6)-(3.8) together with the heat conduction equation, 

kV 2 9- (AT/d)(ve), 



89 
dt 



(3.9) 



and the condition of incompressibility, divv = 0, form the complete set of equa- 
tions which determine the onset of convection in a magnetized ferrofluid. 

Pass over to dimensionless variables by choosing a unit of length d, time 
d 2 /n, velocity n/d, temperature AT, concentration (AT/T)(j), and field potential 
4ttM(AT '/T)d. Then the above equations take the form 



Pr- 1 ^- = -Vp + V 2 v + Rm G 
at 



Q - I # + ^- j 6 + aD<P 



(3.10a) 
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r)9 

|- = V 2 6>- (ve), divv = 0, (3.10b) 
r)0 f)f) f)<T> 

— = LeV 2 Q + ^_ _ ((T - A ) D— + # H (ve), (3.10c) 




fiV^+ia- n)'D$= -D fl- yP+^j 9 > V^ c = 0. (3.10d,e) 

Here Pr = ry/pK is the Prandtl number, Le = D/k the Lewis number, D = d/ dz, 

Rm= m T* ■ G = — ^ + J' (3 - lla ' b) 

where i?m we shall call the magnetic Rayleigh number, $ e is the potential of 
field perturbations outside the fluid layer, and Q = ip + \P6 — (a — fi) D<P is 
introduced as an independent variable instead of the concentration perturbation 
ip. It is worth to note that VJ? has the physical meaning of the matter flux 
perturbation. 

Studying the onset of convection in the Benard configuration, one should put 
all perturbations be dependent on time and horizontal coordinates as exp[At + 
t(k x x + k y y)]. Then from Eqs. (3.10) it follows the system of equations for z- 
dependent amplitudes of z-component of velocity W = (ve), temperature 9, 
matter flux potential J?, and magnetic field potentials and <P C : 

(D 2 -k 2 ){D 2 -k 2 -\/Pr)W-k 2 RmG[f2-(>P + (,£'/£)9 + <TD<P} = , (3.12a) 

(D 2 - k 2 - X)9 -W = , (3.12b) 
ie(D 2 - k 2 - X/Le)Q + MO -(a- fi)\D<P + >P H W = , (3.12c) 
(crD 2 -fik 2 )<P= -D[f2- (& + ££'/£)9], (D 2 - k 2 )<P c = . (3.12d,c) 

4 Boundary Conditions 

To completely specify solutions of Eqs. (3.12), one needs ten boundary conditions 
at confined surfaces with the addition of two conditions on # e far from the layer. 

For the case of rigid boundary surfaces whose heat conductivity is very high 
in comparison with the enclosed ferrofluid, the boundary conditions on velocity 
and temperature read: 

W = DW = 9 = at z = ±1/2. (4.1) 

It is natural to assume such boundaries to be at the same time completely imper- 
vious. The latter condition is expressed by setting the matter flux perturbation 
equal to zero, 

D^ = at 2 = ±1/2. (4.2) 

Boundary conditions on the magnetic potential <P are pretty complicated by the 
fact that <P induces outside the layer a magnetic potential # e governed by Eq. 
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(3.12e). This equation have the solution # e oc exp(±fcz) which must decay far 
from the layer and be valid everywhere outside the layer including the layer 
boundaries too. Hence we come to the relations 

D<P C = -k$ e at z = 1/2 ; D<P C = k$ e at z = -1/2 . (4.3) 

Boundary conditions of continuity of the tangential components of magnetic field 
and the normal component of magnetic induction on the layer boundaries have 
the form 

d<P/dx = d$ e /dx , d$/dy = d<P e /dy , (4.4) 

and 

/ dM dM dM \ 
D * + te (^ + ^ + ^ D *) =D# *' <"- 5 > 
The conditions (4.4) are satisfied if 

<P = $ C at z = ±1/2, (4.6) 

whereas Eq. (4.5) on the substitution M = M(</>, T, £) takes the dimensionless 
form 

aB<P + Q- (w+^j9 = B<P c at z = ±1/2 . (4.7) 

Relations (4.3) and (4.6) permit to write down: D<£ e = =f/c<£ at z — ±1/2, 
that gives us an opportunity to eliminate D^ c from (4.7). Taking into account 
the boundary conditions on the temperature (4.1), we obtain, finally, the closed 
boundary conditions on <P: 

erD<£ + fc<£= -Q at = 1/2; aD<P - k$ = - Q at z=-l/2. (4.8) 

Thus, there is no necessity at all to find the magnetic potential outside the layer. 
The relations (4.1), (4.2), and (4.8) represent the complete set of the rather 
realistic boundary conditions for the set of equations (3.12a)-(3.12d). 
Below, however, we consider the case of idealized boundary conditions 

W = D 2 W = 9 = Q = B<P = at z = ±l/2, (4.9) 

which would be true for free, perfectly heat conducting and pervious boundaries 
made of a superconductor. In the latter case, magnetic field perturbations do 
not penetrate into the bulk of constrained horizontal plates, so that D<£ c = 0. 
Thus, as it follows from (4.7) with allowance for the boundary conditions (4.9) 
on 6 and Q, the value at z — ±1/2 appears to be zero too. 



5 Exact Solution for Free Boundaries 

In the case of boundary conditions (4.9) there exists an analytical solution of 
Eqs. (3.12a) — (3.12d). Satisfying the boundary conditions and the equations by 
means of 

W oc oc fi oc cos 7rz , (p oc sin irz , 
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(tt 2 + k 2 + A)(tt 2 + k 2 + \/Pr){ir 2 + k 2 )[\(air 2 + k 2 ) + Le{f3n 2 + k 2 ){ir 2 + k 2 )} 

-k^RmG \\{£C'/C) + (tt 2 + k 2 + A)<F H + Le{Tr 2 + k 2 ){W + £,£'/£)] = , (5.1) 

where a = jl/ji and f3 = cr/ji. For stationary convective instability (A = 0) the 
magnetic Rayleigh number is given by 

Le {f3ir 2 + k 2 )(ir 2 + k 2 ) 3 



Rm s 



k 4 G [!Ph + Le + iC'/C)} ' 



As it mentioned above, magnetic colloids are characterized by extremely small 
Lewis numbers (typically Le <~ 1CP 4 ). Therefore, using the definitions (2.10) and 
(3.11b), one may rewrite the last equation in the form 

Rm- = Lea 2 {p, 2 + k 2 )j, 2 + k 2 f 

Ak^C 2 {9 + - (A - i)£C/£] ' 

The magnetic Rayleigh number according to its definition (3.12a) is always pos- 
itive. Therefore, as it seen from (5.2), stationary instability sets in only outside 
the interval 

i.e., cither for a large enough positive separation ratio or for a sufficiently small 
negative At the same time, as it will be shown below, inside the interval it 
can occur only oscillatory convective instability. In particular, the oscillatory 
instability is only possible in the absence of the Soret effect, i.e., when 9 = 0. 

Omitting in (5.1) small terms with the coefficient Le, we obtain the cubic 
equation for A: 

(an 2 + k 2 )(ir 2 + fc 2 )A 3 + (Pr + l)(mr 2 + fc 2 )(^ 2 + fc 2 ) 2 A 2 

+ [(mr 2 + fc 2 )(7r 2 + k 2 f - k A Rm G (% + ££'/£)}Pr\ 

= fc 4 (7r 2 + k 2 )PrRmG<P H . (5.4) 

Substituting in this equation A = iw we arrive at two relations which determine 
the magnetic Rayleigh number for oscillatory instability, 

= (Pr + l)a 2 (an 2 + k 2 )(n 2 +k 2 r 

4irflPrk4£ 2 (V + ££>/£)[$ + {££> /£){! + a/fiPr)] ' 1 ' ' 

and the frequency of neutral convective oscillations, 

2 ^ [^-(A-i)6C/AK^ 2 + fc 2 ) 2 

W +{^£' /£){! + a/ fiPr) ' K ™> 

Since Rm osc and uj 2 must be both positive, the oscillatory instability can occur 
only if ^ lies within the interval (5.3). The corridor of oscillatory instability for 
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a ferrofluid with the initial magnetic permeability /io = 3 is plotted in Fig. 1 
for a wide enough range of the Langevin parameter £. On the upper border of 
the corridor, the frequency of convective oscillations turns into zero, whereas on 
its lower border the critical Rayleigh number turns into infinity. In the absence 
of magnetophoresis the corridor grows narrow so that its upper border passes 
along the line \P = 0. In this case the oscillatory instability would be possible 
owing to the Soret effect only for IP < 0. At last, if neither magnetophoresis nor 
Soret effect are in operation (i.e., no particle diffusion), then the cubic equation 
(5.4) for A is reduced to the square one, 

A (an 2 + k 2 ) (tt 2 + k 2 ) [A 2 + (Pr + 1) (tt 2 + k 2 ) A + Pr (tt 2 + k 2 ) 2 ] 

= ArrPrk 4 Rm(^C') 2 . (5.7) 

This equation has no solutions of kind of A = iui. Thus, if there is no diffusion, 
then there is no oscillatory instability as well. For stationary instability Eq. (5.7) 
gives 




2 4 6 8 10 

Undimensional field % 



Fig. 1. The corridor of oscillatory instability (dashed) in (£ ,&) plane. Straight lines 
correspond to \P = 0, & = —0.6, x t = +0.6 and £ = 2. Cross sections of the surface 
Rm c (£,W) along these lines are presented below in Figs. 2-5 respectively. 

Below we shall consider Rm as a function of £ and Any pair of fixed mag- 
nitudes of the parameters determines a neutral stability curve in the (Rm , k) 
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plane. The minima of these curves, Rm c , and the critical circular frequency of 
convective oscillations, / = ui c /2ir, are presented in Figs. 2-5. Figures 2-4 show 
Rm c and / as functions of £ for three magnitudes of & marked in Fig. 1 by 
horizontal lines & = 0, \P = —0.6, and <F = +0.6. Figure 5 represents a stability 
diagram in the plane (Rm c , \P) for the magnitude £ = 2 marked in Fig. 1 by the 
vertical line. 



6 Discussion of Results and Conclusion 

At first, let us consider the case <P = 0. Then Eqs. (5.5)-(5.6) become somewhat 
simpler: 

(Pr + 1) a 2 (an 2 + k 2 ) (ir 2 + k 2 f 



Rm osc = 



Ank A (££') 2 (Pr fi + a) 



2 = 3Pr(y -l)£ 2 (n 2 + k 2 ) 2 

Prfi + a ' 1 ' ' 

The function Rm osc (k) has a minimum Rm° sc at the critical wave number k c 
determined by the formula 

k 2 = (tt/2) 2 (1 - a + ^/l + Ua + a 2 ) . (6.2) 

Here a = jl/fj, < 1 depends on the dimensionless magnetic field £ and the initial 
magnetic permeability /xq. For commonly used ferrofluids (with 1 < ^ < 5) k c 
weakly depends on £. For instance, at fiQ = 3 the critical wave number lies inside 
the interval 2.97 < k c < 3.14 for any £, so that fc c variations do not exceed 6%. 

Substituting k 2 from (6.2) into (6.1) we find the critical Rayleigh number 
and the frequency of neutral oscillations. The values log Rm° sc and / = ui c /2ir 
as functions of £ are presented in Fig. 2. In the same figure it is shown the 
logarithm of Rmi°^ determined by Eq. (5.8) under the replacement k 2 by k 2 
from (6.2). 

We dare say that both the critical Rayleigh numbers, Rm° sc and Rmi°^ 
have a direct physical meaning. Depending on conditions of a real experiment 
it will set in either stationary convection above Rm^ or oscillatory convection 
above Rm° sc . The point is that ferrofluids are characterized by two very differ- 
ent time scales [13]. One of two is the above-mentioned long mass diffusion time 
to = d 2 /ir 2 D: even in a thin layer (~ 1 mm) the concentration equilibrium is es- 
tablished for a time longer than one hour. Another one is a short thermodiffusion 
time tt = d 2 /tt 2 k: in the same fluid layer the stationary temperature gradient 
sets in merely in one second. The ratio tt/td is equal just to the Lewis number 
Le = D/k ~ 10~ 4 . The presence of such different time scales makes possible two 
scenarios of the onset of convection. First, the temperature difference between 
the layer boundaries can be increased from zero till a supercritical value during 
a time r <C td- Then diffusion processes have no time to evolve, so the magnetic 
colloid behaves like a pure fluid. Indeed, magnetic grains do not move without 
diffusion relative to the surrounding liquid matrix, they are "frozen" into it. Just 
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T 




Fig. 2. Critical magnetic Rayleigh number Rm c and circular frequency of neutral os- 
cillations / versus magnetic field strength £. Curve (a) - oscillatory solution Rm° sc 
obtained with allowance for magnetophoresis but without Soret effect = 0), curve 
(b) - stationary solution Rm^ in the absence of the particle diffusion (i.e., there are 
neither magnetophoresis nor the Soret effect). 



that very case has been studied in [3-8] . Stationary convection starts in this case 
at Rmi°^ . The latter value has a minimum at a certain £ ; for a water-based 
ferrofluid (Pr = 7) with fi a = 3 this minimum, min[i?mc°^] = 1742, is reached 
at £ = 2.52 and k c = 3.02. 

Second scenario can be realized in a thin enough fluid layer, where the long 
diffusion time, td, is not too long. Then it is possible to adjust so low heating 
rate, that the temperature difference across the layer, AT, will be formed for 
a time r > td- Under this condition the concentration gradient induced by 
the temperature gradient due to magnetophoresis (the Soret effect is absent at 
& = 0) will be built up undisturbed by convection. In this case, as it shown above, 
there arises oscillatory instability at Rm° sc . For a ferrofluid under consideration 
the critical Rayleigh number reaches the minimum, min [Rm° sc ] = 4993, at £ = 
1.47; corresponding wave number and frequency of neutral oscillations are k c = 
3.08 and u> c = 12.13. The frequency increases monotonously with £ approaching 
the saturation value w c = 34.19 at £ — > oo. 

Let us estimate the temperature difference and other parameters required 
to cause oscillatory convection by the thermomagnetic mechanism. Assuming 
the mean diameter of magnetite particles (M s = 480 G) to be 2a = 11 nm, we 
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obtain the particle magnetic moment m = (Ana 3 /3)M S = 3.345 x lCT 16 erg/Oe, 
so the value /io = 1 + 4ir(f>M s m/3kBT = 3 is reached at the room temperature 
at <f> = 12%, and £ = mH/k B T = 1 at H = 122 Oe. Thus the minimum in 
the curve of log Rm° sc in Fig. 2 takes place at H = 180 Oe (to get such a field 
within the ferrofluid layer one must apply externally H e = /J-oH = 540 Oe). The 
minimum Rayleigh number (« 5000) is reached at (AT) C = (1.02/e?) K/cm, thus 
the critical temperature difference is about 10 K in the layer of d = 1mm. For this 
layer the unit of time, cP/k, is equal to 6.9 s, so the above-pointed dimensionless 
frequency ui c = 12.13 of critical convective oscillations corresponds to the period 
of 3.6 s. 

Maps of stability in the (Rm c , £)-plane are presented in Figs. 3 and 4 for 
<P = —0.6 and ^ = +0.6, respectively. As it seen in Fig. 1, for both these 
magnitudes of the separation ratio there are certain intervals of £ in which it can 
occur oscillatory instability. Out of the intervals there arises stationary Soret 
convection which is characterized by very small Rayleigh numbers (5.2) due to 
the smallness of the Lewis number (for a water-based ferrofluid Le = 1.9 x 10~ 4 ). 
The critical Rayleigh number for stationary instability, Rmf", is obtained from 
(5.2) under the substitution k% from (6.2) with the replacement of a by fi. 




Fig. 3. Critical magnetic Rayleigh number and frequency of neutral oscillations versus 
magnetic field strength. Curve (a) - oscillatory solution Rm° BC and (b) - stationary 
solution Rmf with allowance for magnetophoresis and the negative Soret effect (<P = 
—0.6); curve (c) - stationary solution Rmi°^ in the absence of the particle diffusion. 
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Fig. 4. Critical magnetic Rayleigh number and frequency of neutral oscillations versus 
magnetic field strength. Curve (a) - oscillatory solution Rm° sc and (b), (c) - stationary 
solution Rmf with allowance for magnetophoresis and the positive Soret effect (\P = 
0.6); curve (d) - stationary solution Rm£^ in the absence of the particle diffusion. 



For a fixed negative separation ratio, > W > — 1, oscillatory instability 
takes place within an interval of the field < £ < £o (see Figs. 1 and 3), where 
£o satisfies the equation 

# + 6C'(0/£ = 0. (6.3) 

On the borders of the interval, the critical Rayleigh number (5.5) turns into 
infinity as for £ — > and as (£ — £o) _1 for £ ^ £ (see Fig. 3). The frequency 
of neutral oscillations, w c (£), commences from a finite magnitude at £ = 0, 



w c (0) = 2tt 2 



Pr |tf| 



1 + Pr(l- \&\) 



and grows monotonously with £ until it reaches £o- 

The upper border of the interval (5.3) is described by 

M 

where fi = l + 3(yU — !)£'(£)• The function (6.4) has a maximum $ m at a certain 
value £ m dependent on the initial magnetic permeability /i = 1 + 47rx - For the 
fluid under consideration (^ = 3) we find £ m = 3.90 , & m = 0.806. For any 
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positive \P < >P m there exists an interval £1 < £; < £2 of oscillatory instability - 
see Figs. 1 and 4. The frequency of oscillations reaches a maximum within the 
interval and turns into zero on its borders. Oscillatory branch (a) of iim c (£) links 
in Fig. 4 two codimension-2 points, i?m c (£i) and Rm c (^2), in which stationary 
branches (b and c) bifurcate. 

A diagram of stability in the (Rm c , i^-plane at the fixed magnitude of the 
field, £ = 2, is shown in Fig. 5. This diagram looks more or less traditional for 
binary mixtures [14,15] with exception of two things. Firstly, the codimension-2 
point (i.e., the intersection point of a and c branches) is usually located in the 
close vicinity of & = 0. In Fig. 5 this point has essentially shifted towards positive 
\P due to magnetophoresis. Secondly, the stationary branch & as a general rule is 
located in the region of negative Rayleigh numbers, i.e., the branch determines 
the onset of convection in the binary mixture heated from above. Thermomag- 
netic mechanism of convection knows, however, neither the top nor the bottom: 
the magnetic Rayleigh number (3.11a) is proportional to (AT) 2 , so it is always 
positive. 




Fig. 5. Critical magnetic Rayleigh number and frequency of neutral oscillations versus 
separation ratio at the fixed magnetic field strength £ = 2. Curve (a) - oscillatory 
instability, curves (6) and (c) - stationary instability. 



Thus, convective instability of magnetized ferrofluids is strongly effected by 
the magnctophoretic and thermophoretic (the Soret effect) transfer of magnetic 
grains. The most interesting consequence of the transfer is an opportunity to ob- 
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Fig. 6. Stationary instability curves for different /lio in the case when the particle 
diffusion is not operative. 

serve oscillatory convective instability. However, as we have already commented, 
to observe predicted oscillations one needs to satisfy certain experimental con- 
ditions. Namely, the temperature difference must not be increased faster than 
the limit imposed by concentration diffusion. In the opposite fcrrofluid 
behaves like a pure fluid what results in stationary instability at the critical 
Rayleigh number Rm^ ; its dependence on £ for some magnitudes of is shown 
in Fig. 6. In conclusion, it is worth to remind that all our results are relevant 
to thin fcrrofluid layers (~ 1mm) where magnetic mechanism of convection pre- 
dominates over the buoyancy mechanism, and the characteristic diffusion time, 
td, is not too long. To compare above-stated theoretical predictions with future 
experimental results one should perform a more precise analysis, which would 
satisfy realistic boundary conditions (4.1), (4.2) and (4.8) instead of idealized 
boundary conditions (4.9) used in the paper. 
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Glossary 

(f> concentration of magnetic grains (volume fraction) 

M magnetization 

H magnetic field 

B magnetic induction: B = H + 47rM 

F magnetic force (volume density) 

v ferrofiuid velocity 

T temperature 

fee Boltzmann's constant 

p mass density 

rj fluid viscosity 

Xo initial magnetic susceptibility 

fiQ initial magnetic permeability: fiQ = 1 + ^Xo 

X magnetic susceptibility: x = M/H 

fi magnetic permeability: /j, = 1 + 4ir(M/H) 

D mass diffusion coefficient 

k thermal diffusivity 

5t Soret coefficient 

& separation ratio for magnetic convection: 9 = (T/0)St 

t time 

td characteristic diffusion time 

tt characteristic thermodiffusion time 

Le Lewis number: Le — tt/td 

Pr Prandtl number 

Ra Raylcigh number 

Rm magnetic Rayleigh number 

k wave number 

£ Langevin's parameter (dimensionless magnetic field) 

£(£) Langevin's function: £(£) = coth^ — 

A increment of a normal mode 

uj angular frequency of neutral oscillations 

/ frequency of neutral oscillations: / = u>/2tt 



